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Abstract 

We have designed and carried out simulation studies of a two layer Shower 
Maximum Detector diagonally off-set (SMD-dos) optimized for the separa- 
tion of 7T° showers from 7 showers in the 30 to 150 GeV energy range. For 



70, 



90% 7 acceptance the SMD-dos yields ir rejection efficiencies of 92±4 
87±4 % and 32±2 %, respectively, for 30, 50 and 150 GeV incident energies. 
We find that the SMD-dos is superior to a conventional geometry single-layer 
or mutiple-layer shower maximum detector (SMD), of equal granularity, by 
an average factor of ~ 1.5 over the 50 to 150 GeV energy range. We also 
find, the unexpected result, that the SMD-dos gives better 7r° rejection, for 
the same number of channels, than a SMD. At hadron - hadron colliders the 
signature of choice for the detection of the Higgs particle, in the mass range 
of 120 to 160 GeV, is via the decay H — > 77. The addition of a SMD-dos to 
the planned detectors at the LHC would significantly reduce the background 
to the 7 signal coming from prolific 7r° production. 



1 Introduction 



The challenging task of separating ir° initiated showers from 7 initiated show- 
ers, at high incident energies, has over the past few years attracted a lot of 
interest and motivated extensive R and D work |], 0, [3|, f|, ||, [7j . The sched- 
uled building of the LHC machine has added new interest to the ■y/n sepa- 
ration problem. The main motivation for building the LHC machine is the 
quest for the Higgs particle. The Standard Model has weathered extremely 
well strenuous experimental scrutiny during the past decades; however, it still 
lacks experimental support for the predicted mass generating Higgs mecha- 
nism. The only viable explanation of the masses of the elementary particles in 
the framework of the Standard Model is the Higgs mechanism. The two main 
collaborations approved for operation at the LHC, ATLAS || and CMS|| 
state in their Technical Reports the paramount importance of the searches 
for the Higgs particle at the LHC machine. The upgrading of LEP to LEP2 
will allow the search for the Higgs particle up to a mass of about 100 GeV. 
Searches for the Higgs particle at higher masses will have to wait for the 
operation of the LHC machine. 

We have designed and carried out simulation studies of a two layer Shower 
Maximum Detector diagonally off-set (SMD-dos) optimized for the separa- 
tion of 7T° showers from 7 showers in the 30 to 150 GeV energy range. A 
comparison of the SMD-dos with a conventional geometry single-layer and 
mutiple-layer shower maximum detector (SMD), having the same granular- 
ity, shows that the SMD-dos yields higher ir° rejection efficiency in the energy 
range studied. 

We have also carried out a comparison of the SMD-dos versus a single 
layer SMD with both detectors having the same number of channels. We find 
the unexpected result that the SMD-dos yields better 7r° rejection efficiency 
than the SMD. 

The electromagnetic showers were generated and studied using the Monte 
Carlo program GEANT 3.15 implementing a 100 KeV cut for 7's and elec- 
trons. The shower maximum detector, either a SMD-dos or a SMD, was 
embedded inside a E.M. calorimeter, consisting of a mixture of Pb and scin- 
tillator (CH) having a radiation length of 0.8 cm. The performance of the 
SMD-dos was optimized by carrying out extensive Monte Carlo studies of 
the 7/71"° separation efficiency as a function of the following design parame- 
ters: number of scintillator layers, thickness of scintillator layer, granularity 
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(cell size) and position of the SMD-dos or SMD inside the E.M. calorimeter. 
Typically, samples of 500 showers were generated for each configuration of 
design parameters and incident energy. The distance of the E.M. calorimeter 
from the interaction region was kept fixed at 150 cm, the design distance of 
the CMS detector approved to operate at the LHC machine, and the showers 
were generated with the shower initiating particle, a 7 or 7r°, normally inci- 
dent on the E.M. calorimeter. We have studied the stand-alone capabilities 
of the SMD-dos; however, since the algorithms developed assume that the 
energy of the shower, in question, is known, the SMD-dos must be used with 
an E.M. calorimeter to determine the total energy of the shower. 

At hadron - hadron colliders the signature of choice for the detection of 
the Higgs particle, in the mass range of 120 to 160 GeV, is via the decay 
H — > 77 [[Kj]. These searches will encounter high 7 backgrounds coming 



from prolific tt° production which decay into two 7's. The attainment of 
good 7/71- separation is of crucial importance to the success of the proposed 
searches. The incorporation of a SMD-dos into the planned detectors at the 
LHC would appreciably reduce the 7 background. 



2 The Geometry of the Diagonally Off-set 
Shower Maximum Detector 

The geometry of the SMD-dos, shown in Figure p], consists of 2 identical 
scintillator layers, placed one behind the other. Each layer is subdivided 
into square cells which are individually read out and the energy deposited 
in each cell is recorded. The 2 layers are diagonally off-set with respect to 
each other by half a cell, so that the intersection of any four adjacent cells 
in one layer corresponds to the center of a cell in the other layer. Studies of 
the 7/71" separation efficiency of granular detectors show that the separation 
efficiency decreases as a function of the distance of the shower center from the 
nearest intersection. Hence increasing the number of intersections per unit 
area enhances the separation efficiency. The area of each cell, in the front 
layer, overlaps symmetrically one quarter of the area of each of four adjacent 
cells in the back layer. And likewise the area of each cell, in the back layer, 
overlaps symmetrically one quarter of the area of each of four adjacent cells 
in the front layer. Therefore, the energy collected by a particular cell is 
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subdivided unto 4 cells giving a substantial increase in granularity. 
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Figure 1: Layout of the two layers of the SMD-dos. The cells of the front 
layer are shown in solid lines and the cells of the back layer are shown in 
dashed lines. One back layer cell is highlighted to show the four quarter cells 
in the front layer which overlap it. 

In the singular case that all the energy deposited by a shower is confined, 
in a particular layer, to one cell it is impossible with a conventional one-layer 
or multiple-layer SMD to determine whether the shower was initiated by a 
7 or 7T°. However, the unique geometry of the SMD-dos, which subdivides 
the energy collected by a particular cell, in one layer, unto four cells in the 
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other layer, does furnish some degree of 7/71" separation efficiency even in 
this singular case. 



3 Definition of Variables Used in the Algo- 
rithms 

The energy and position information collected, per event, by the scintillator 
cells is used to construct 10 variables. 
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Figure 2: Display of cell configurations used in the determination of E max , 
E4, E 8 and £12. 

These variables embody the experimental information which goes into the 
7/71-° separation algorithms. The information collected in the first layer is 
used to construct 5 variables and likewise the information from the second 
layer is used to construct the other 5 variables. The detailed description and 
construction of these 10 variables is detailed in this section. In Figure are 
displayed the cell configurations used for the determination of E max , £4, Eg 
and E12. 
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• < r >, the average energy weighted radius of the shower. The deter- 
mination of < r > from the data collected by the cells in the layer is 
discussed in the following paragraph. 

• E max , the maximum energy deposited in a single cell, in the layer. See 
Fig. 2a. 

• E 8 , the sum of energies deposited in the eight cells which surround the 
cell with the maximum energy. E 8 is normalized to E tota i. See Fig. 2b. 

• £4, the sum of energies deposited in a 4 cell square which includes the 
cells with the maximum and next to maximum energy. £4 is normalized 
to E total . See Fig. 2c. 

• Ei2i the sum of energies deposited in the 12 cells which surround the 4 
cell square described in the preceding item. E\ 2 is normalized to E tota i ■ 
See Fig. 2d. 

The energy information gathered by each cell is read out separately and 
is labelled, Eij, where the x and y coordinates of the particular cell, in 
the two dimensional layer, are denoted, respectively, by the i and j sub- 
scripts. The average energy weighted radius, < r >, is calculated for each 
layer separately. First the coordinates of the energy weighted center of 
the shower, xo and yo, are calculated as follows: xo = J2i,j E^Xij J2ij E^ 
and yo = EijVi/ J2i,j Eij, where the i and j summations are carried out 
over all the cells that have collected energy and for i < j, to avoid dou- 
ble counting. Using the x and y values, we determine the distances 
between the center of the ij cell and the shower energy weighted center, 



r ij = y ( x o ~ x i) 2 + (Vo ~ Vj) 2 - T ne avera g e energy weighted radius of the 
shower, < r >, is determined as follows: 



Information from both layers is used to construct the 2 layer correlation 
variables. In the following we differentiate between the data gathered in the 
layer that recorded the cell with the maximum energy, the leading layer, and 
label them with the superscript, /, and the data from the other layer, the 
non-leading layer, and label them with the superscript, n. 




<r >= 
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The following 6 variables embody the 2 layer correlation information used 
in the algorithms: the correlations between E max from one layer and £4 from 
the other layer are contained in the variables C\ and C 2 . Similarly the 
correlations between E 4 and E 8 are contained in the variables C% and C4 
and, finally the correlations between E 8 and E 12 are contained in C 5 and C 6 . 

• C\ is the difference, E2 — E l max , normalized to E to tai- 

• C2 is the difference, E\ — E^ ax) normalized to E tota i- Notice that C2 is 
similar to C\ with the roles of the 2 layers interchanged. 

• C 3 is the difference E% — E l 8 normalized to E tota i- 

• C 4 is the difference E\ — Eg normalized to E tota i. C 4 is similar to C 3 
with the roles of the 2 layers interchanged. 

• C 5 is the difference Eg — E[ 2 normalized to E tota i. 

• Cq is the difference E\ — E™ 2 normalized to E to tai- C% is similar to C5 
with the roles of the 2 layers interchanged. 

4 Analysis 

The study of the separation of 7/71"° showers is quantified by determining the 
rejection efficiency of tt° initiated showers as a function of the acceptance 
of 7 initiated showers. Eight variables are used in these analyses: the six 
correlations, Ci,.. ,Cq and the two average energy weighted radii, < r > l and 

< r > n . 

For every configuration of design parameters and incident energy stud- 
ied, two samples of showers, one initiated by 7's and the other by 7r°'s, were 
generated and propagated through the E.M. calorimeter and the SMD-dos. 
The shower information collected by the SMD-dos was used to construct 
the 8 variables previously described. A program which uses the Simulated 
Annealing Algorithm was utilized to delineate an 8-dimensional volume, cor- 
responding to the 8 variables used, which contains a fixed percentage of the 
7 initiated showers while keeping to a minimum the number of n° initiated 
showers. The delineation of this volume sets upper or lower cutoff values for 
each of the 8 variables. The 7r° rejection efficiency is equal to the percentage 
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of 7T° initiated showers which are excluded from the delineated 8-dimensional 
volume. 

The shower information collected by the SMD was also processed using 
the same analysis chain detailed above. However, since for the SMD no 
correlation information is available, the following variables were used: E4, 
E 8 , E 12 and < r >. 

The computer flow chart of the program used along with a short descrip- 
tion of those aspects of the program which are specific to the 7/71° separation 
problem can be found elsewhere @. In Fig. ||] are shown the projections, for 
50 and 100 GeV, of the 8-dimensional volume of 7 initiated showers and 
7T° initiated showers, separately, unto the plane spanned by the correlation 
variables C\ and C3. 
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Figure 3: Scatterplots of the correlation variables C3 versus C± for showers 
initiated by 7's and 7r°'s. a) 50 GeV 7 showers and b) 50 GeV 7T° showers, c) 
100 GeV 7 showers and d) 100 GeV 7r° showers. The respective cutoff values 
of Ci and C 3 , for 90% 7 acceptance, at 50 GeV and 100 GeV, are shown by 
dashed rectangles on the lower right corner of each figure. 

The cutoff values of C\ and C 3 for 90% acceptance of the 7 initiated 



7 



showers define the rectangles shown in the lower right corners. For 50 GeV 
incident energy the 7r° initiated showers overwhelmingly fall outside the ac- 
ceptance rectangle giving high ir° rejection efficiency while for 100 GeV inci- 
dent energy the number of tt° initiated showers excluded from the acceptance 
rectangle is smaller hence giving lower ir° rejection efficiency . 

The cutoff values for different incident energies are given in Table |I[ 
These cutoff values as well as all other results presented in this paper, for 
the SMD-dos and SMD, were determined from samples of 500 showers each 
which were generated under the following conditions: 

• The distance of the E.M. calorimeter from the interaction region was 
fixed at 150 cm. 

• The shower initiating particle, 7 or ir°, was incident normally on the 
E.M. calorimeter. However, the direction of the initiating particle was 
spread randomly within a small angular range corresponding to an area 
of a few cells of the SMD-dos or SMD. 

• The 7T° rejection efficiency was calculated for 90% 7 acceptance. 

• The SMD-dos or SMD was positioned at a depth of 6 X^s inside the 
E.M. calorimeter. 

• Thickness of the scintillator layers set at 4 mm. 

• The nominal granularity, cell size, is 1 x 1 cm 2 . In the event that a 
different cell size was used, it is explicitly stated. 



E- 

J -'in 


Ci< 


C 2 < 


c 3 > 


c 4 > 


c 5 < 


c 6 < 


(r) 1 < 


(r) n < 


30 


0.18 


0.70 


0.65 


0.21 


0.65 


0.17 


0.58 


0.80 


50 


0.17 


0.65 


0.67 


0.22 


0.63 


0.15 


0.55 


0.82 


80 


0.17 


0.68 


0.70 


0.20 


0.62 


0.16 


0.53 


0.85 


100 


0.16 


0.70 


0.71 


0.22 


0.65 


0.18 


0.50 


0.88 


120 


0.15 


0.67 


0.71 


0.23 


0.60 


0.14 


0.50 


0.90 


150 


0.15 


0.64 


0.73 


0.23 


0.58 


0.15 


0.48 


0.91 



Table 1: Cutoff values of the six correlation variables, Ci,.. ,Cq, and the two 
average energy weighted radii, < r > l and < r > n , for incident energies of 
30 to 150 GeV. 
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5 Results and Discussion 



The 7T° rejection efficiency as a function of the incident energy of the show- 
ering particle is shown in Fig. |]a for both the SMD-dos and SMD. Both 
detectors have equal granularity, lxl cm 2 cells. The 7r° rejection efficiency 
of the SMD-dos and the SMD decrease as a function of energy; however, the 
7T° rejection efficiency of the SMD-dos is higher than that of the SMD over 
the whole energy range. The ratio of the rejection efficiency of the SMD- 
dos to the SMD, the 'advantage factor', is plotted as function of energy in 
Fig. |b. 
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Figure 4: a)The ir° rejection efficiency of the SMD-dos (•) and a conventional 
SMD (a) as a function of incident energy. b)The 'advantage factor' (■) of 
the SMD-dos over the SMD as a function of incident energy. 



9 



The value of the 'advantage factor' is close to 1.0 at 30 GeV and it then 
rises to a maximum value of 1.63 at 80 GeV. With increasing energy the 
value of the 'advantage factor' decreases gradually reaching the value of 1.39 
at 150 GeV, the highest energy studied. In Table || are summarized the 
the 7T° rejection efficiencies of the SMD-dos and SMD, of equal granularity, 
for various incident energies. The 'advantage factor' has been calculated for 
each incident energy. The following parameters, which are relevant to 7/71-° 
separation problems are given, per energy, in Table ^|: 612, the minimum 
opening angle in the lab system of the n° — > 77 decay; c?i j2 , the distance 
between the two 7 centers at the position of the SMD-dos or SMD and the 
ratio di^/A, where A is the length of cell side, 1 cm. 



E 

'-'in 


7T° Rejection 


e 12 


di,2 


di,2 


7T° Rejection 


'advantage 


[GeV] 


SMD-dos [%} 


[mrad] 


[cm] 


A 


SMD {%] 


factor' 


30 


92±4 


9.3 


1.43 


1.43 


87±4 


1.06±0.07 


40 


90±4 


7.0 


1.07 


1.07 


80±4 


1.13±0.08 


50 


87±4 


5.6 


0.86 


0.86 


62±3 


1.40±0.09 


80 


65±3 


3.5 


0.55 


0.55 


40±2 


1.63±0.11 


100 


54±3 


2.8 


0.43 


0.43 


34±2 


1.59±0.13 


120 


39±2 


2.3 


0.36 


0.36 


26±2 


1.50±0.14 


150 


32±2 


1.9 


0.29 


0.29 


23±2 


1.39±0.14 



Table 2: The 7r° rejection efficiency of the SMD-dos and SMD, of equal 
granularity, is presented for different values of incident energy. The following 
parameters, which are explained in the text, are given per incident energy: 
©12, ^i,2) the ratio cii^/A and the 'advantage factor'. 



The dependence of the ir° rejection efficiency of the SMD-dos and SMD 
on the granularity of the detector was studied at 50 and 100 GeV incident 
energies. The results of studies are shown in Figure |5] and are summarized 
in Table |3|. The expected increase in ir° rejection efficiency with increasing 
detector granularity for both SMD-dos and SMD is seen in Fig. 5a. 

For 50 GeV incident energy the 7r° rejection efficiency of the SMD-dos 
and the SMD converge as the granularity increases. The 'advantage factor', 
plotted in Fig. 5b, is seen to fall rapidly from a value of 1.71 for granularity 
of 7 x 7 to a value of 1.09 for a granularity of 14 x 14. The 'advantage factor' 
for granularity of 20 x 20 is consistent with 1.0. For 100 GeV incident energy 
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the 7T° rejection efficiency of both the SMD-dos and the SMD increase with 
increasing granularity; however, the 7r° rejection efficiency of the SMD-dos 
is higher over the whole energy range. The 'advantage factor' has a value of 
1.4 for 7x7 granularity, it attains the maximum value of 1.7 for granularity 
14 x 14 and it then decreases to 1.5 for 20 x 20 granularity. 
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Figure 5: a) The 7r° rejection efficiency, as a function of detector granularity, 
of the SMD-dos (a) and a conventional SMD (□) at 50 GeV incident energy; 
at 100 GeV, SMD-dos (a) and a conventional SMD (■). b) The 'advantage 
factor' as a function of detector granularity at 50 GeV (o) and 100 GeV (•). 
The granularity is that of a 10 x 10 cm 2 layer divided into n x n square cells. 
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From the data given in Table |3| it is possible to compare the ir° rejection 
efficiency performance of the SMD-dos versus the SMD for the same number 
of channels. For 50 GeV incident energy a SMD-dos detector which has 
two 7x7 layers needs 98 channels in order to read out its 98 cells while a 
SMD detector of 10 x 10 granularity needs 100 channels. From Table we 
determine that the SMD-dos is a factor of 1.16 better, 72% vs 62%. A 10 x 10 
SMD-dos is better than a 14 x 14 SMD by a factor of 1.07. For 100 GeV 
incident energy we find that a 10 x 10 SMD-dos is superior to a 14 x 14 SMD 
by a factor of 1.38. Also a 14 x 14 SMD-dos is a factor of 1.18 better than a 
20 x 20 SMD. 



Energy = 50GeV 


Number 


A 




7r°Rejection 


7T° Rejection 


'advantage 


of cells 


[cm] 


A 


SMD-dos [%] 


SMD [%] 


factor' 


7x7 


1.43 


0.60 


72±3 


42±2 


1.71±0.11 


10x10 


1.0 


0.86 


87±4 


62±3 


1.40±0.09 


14x14 


0.72 


1.20 


88±4 


81±4 


1.09±0.07 


20x20 


0.50 


1.72 


91±4 


88±4 


1.03±0.07 


Energy = 100 GeV 


Number 


A 




7T° Rejection 


7T° Rejection 


'advantage 


of cells 


[cm] 


A 


SMD-dos [%] 


SMD [%] 


factor' 


7x7 


1.43 


0.30 


36±2 


26±2 


1.38±0.13 


10x10 


1.0 


0.43 


54±2 


34±3 


1.59±0.13 


14x14 


0.72 


0.60 


67±2 


39±3 


1.71±0.12 


20x20 


0.5 


0.86 


87±3 


57±4 


1.52±0.11 



Table 3: The 7r° rejection efficiencies of the SMD-dos and SMD, the 'ad- 
vantage factor's, the values of A and the ratio di^/A, are given for different 
detector granularities, for 50 and 100 GeV incident energies. 

The advantage of using 2 layers which are diagonally off-set with respect 
to each other as oppossed to using a single layer, given that the total number 
of cells of each detector is the same, can be understood as follows: the j/n 
separation power of a granular detector is strongly dependent on the distance 
of the shower center from an intersection. The success rate of any available 
7/71-° separation algorithm is higher when the deposited shower energy is 
distributed among many cells. As previously stated, if all the deposited 
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energy is collected by a single cell it is impossible to determine the identity 
of the showering particle. The distance between intersections for a SMD-dos 
is smaller by a factor of 1/ \[2 than that of a SMD when both detectors have 
equal number of cells. Therefore, a showering particle incident on a SMD-dos 
detector will, on the average, have its shower center nearer to an intersection 
than if it were incident on a SMD detector. 




Figure 6: a) The 7r° rejection efficiency of the SMD-dos (•) and a conventional 
SMD (o) as a function of the ratio, d lj2 /A. b) The 'advantage factor' as a 
function of the ratio, rf 12 /A. 



The ratio of the 2 variables, d± : 2, the distance between the two 7 centers 
at the position of the SMD-dos or SMD and A, the length of cell side, 0^,2/ A, 
is a useful a measure of the 7/71-° separation difficulty The 7/71- separation 



13 



becomes more difficulty as the ratio dip/ A tends to small values while it is 
significantly easier as the ratio tends to values above 1.0. The behaviour of 
the 7T° rejection efficiency and the 'advantage factor' as a function of dip/ A 
are shown in Figure ^. The 7r° rejection efficiency, as expected, increases for 
both the SMD-dos and the SMD as the value of this ratio increases. We 
notice that the 'advantage factor' has an appreciable value only in the range 
of from ~0.2 to ~1. Hence we find that when the 7/71"° separation is very 
difficulty, values smaller than ~0.2, than the efficiencies of the SMD-dos and 
SMD are comparable and likewise when the separation becomes easy, values 
higher than ~1.0, than again no advantage of SMD-dos over SMD. 

If we express, dip in terms of the incident energy, Ei nc , mass of the ir° 
and, L, distance to SMD-dos or SMD from the interaction point, then we 
can write the following expression: 

dip 2 x m w x L 
A E inc x A 

Hence the ratio, dip/ A, is a function of 3 independent variables, A, L 
and Ei nc . So for instance if we fix the value of 2 of these variables we can read 
off from Figure [] the permissible range of the third variable, for a desired tt° 
rejection efficiency. Figure ^| can be used, in general, in order to get the 7r° 
rejection efficiency and 'advantage factor' for different combinations of the 3 
variables not directly studied in this paper. However, it would be prudent to 
limit the use of these figures to values of the 3 variables which are within a 
factor of 2 of the ranges explicitly studied in this paper. 



6 Conclusions 

We have designed and extensively studied a new geometry SMD consisting 
of two layers of scintillator which are diagonally off-set with respect to each 
other. The layers of scintillator are subdivided into square cells. This geome- 
try has the virtue of increasing substantially the number of intersections, per 
unit area, seen by the incoming shower. A comparison of the SMD-dos with a 
conventional geometry single-layer and mutiple-layer SMD, having the same 
granularity, shows that the SMD-dos yields higher ir° rejection efficiencies. 
For 90% 7 acceptance the SMD-dos yields 7r° rejection efficiencies of 92±4 %, 
87±4 % and 32±2 %, respectively, for 30, 50 and 150 GeV incident energies. 
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We find that the SMD-dos is superior to the SMD, of equal granularity, by 
an average factor of ~ 1.5 over the 50 to 150 GeV energy range. We also find 
that the SMD-dos gives better 7r° rejection, for the same number of channels, 
than a SMD. At 100 GeV a 10 x 10 cells per layer SMD-dos, which utilizes 
200 channels, is a factor of ~ 1.4 better than a 14 x 14 cells, single layer SMD, 
which utilizes an equal number of channels. The 7r° rejection efficiency of the 
SMD-dos and SMD and the 'advantage factor' are plotted as a function of 
the ratio, di )2 /A. With the help of these plots it is possible to investigate 
the relative importance of the 3 variables, A, L and E inc , over a range which 
goes beyond the work presented in this paper. 

At hadron - hadron colliders the signature of choice for the detection of 
the Higgs particle, in the mass range of 120 to 160 GeV, is via the decay 
H — > 77. The addition of a SMD-dos to the planned detectors at the LHC 
would significantly reduce the background to the 7 signal coming from prolific 
7T° production. 
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